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Background: Theoretical calculations have shown that the energy and angular correlations in the three-body 
decay of the two-neutron unbound can provide information on the ground-state wavefunction, which has been 
predicted to have a dineutron configuration and 2n halo structure. 

Purpose: To use the experimentally measured three-body correlations to gain insight into the properties of ^®0, 
including the decay mechanism and ground-state resonance energy. 

26 27 24 

Method: "°0 was produced in a one-proton knockout reaction from F and the O +n + n decay products 
were measured using the MoNA-Sweeper setup. The three-body correlations from the ^®0 ground-state resonance 
decay were extracted. The experimental results were compared to Monte Carlo simulations in which the resonance 
energy and decay mechanism were varied. 

Results: The measured three-body correlations were well reproduced by the Monte Carlo simulations but were 
not sensitive to the decay mechanism due to the experimental resolutions. However, the three-body correlations 
were found to be sensitive to the resonance energy of ^®0. A la upper-limit of 53 keV was extracted for the 
ground-state resonance energy of ^®0. 

Conclusions: Future attempts to measure the three-body correlations from the ground-state decay of ^®0 will 
be very challenging due to the need for a precise measurement of the momentum at the reaction point in the 
target. 

PACS numbers: 21.10.Tg, 23.90.-l-w, 25.60.-t, 29.30.Hs 


I. INTRODUCTION 

Exploring the limits of the chart of the nuclides pro¬ 
vides fundamental benchmarks for theoretical calcula¬ 
tions and opportunities to discover new phenomena [II- 
1^. Recently, experimental and theoretical studies of 
two-proton (2p) and two-neutron (2n) unbound systems 
have garnered significant interest [^-Q- These systems 
allow for the properties and decay mechanisms of nu¬ 
clei with extreme neutron-to-proton ratios, existing be¬ 
yond the driplines, to be examined. Of particular in¬ 
terest are the correlations present in the decay of these 
three-body systems which can offer new insights into the 
initial structure and configuration of the nucleus [Ml- 
At the driplines, unique situations can present them¬ 
selves in which the sequential decay process is forbidden 
and the “true” two-nucleon decay can be observed El- 
While substantial progress has been made in the theoret- 
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ical descriptions of the 2p decay mechanism, resulting in 
impressive reproductions of experimental measurements, 
theoretical frameworks for describing the full three-body 
decay of 2n unbound systems are still at an early stage 
and appear to pose new challenges in comparison to the 
2p decay 0, [ffil, [H [H . 


While the correlations from the decays of a wide range 
of 2n unbound systems have been measured (®H M, 

23, i3Li P, 


21 


23), the 260 sys- 


^°He [H, [11, ^^Li (excited state) 

^^Be (excited state) [H, and ^6Be 
tern provides a particularly interesting case to examine 
in light of recent experimental and theoretical work. Cur¬ 
rently, constraints on the ^60 ground-state resonance en¬ 
ergy of < 200 keV and < 40 keV have been reporte d by 
the MoNA collaboration [^l and R3B collaboration [24| . 
respectively. Thus, the sequential decay of ^60 through 
2^0, with a ground state unbound by 770 keV [^[^. 
is forbidden and 260 must decay directly to 2^0 through 
the simultaneous emission of two neutrons. These con¬ 
straints have provided a sensitive observable for ab ini¬ 
tio-type calculations examining the role of three-nucleon 
forces and continuum effects (26l - l^ . Furthermore, the¬ 
oretical calculations of Grigorenko et al. suggested that 
26 O would be a candidate for 2n radioactivity assum¬ 
ing a pure [^ 3 / 2 ]^ neutron configuration and near thresh¬ 
old ground-state resonance [13. Through the use of the 
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decay-in-target technique [3l|, the MoNA collaboration 
extracted a half-life for of 4.5ti'5(stat) ± 3(syst) pi¬ 
coseconds which suggested the possibility for two-neutron 
radioactivity [s^- Based on the calculations of Grig- 
orenko et al. a fi /2 = 4.5 ps w ould only require the decay 
energy to be < 100 keV [^. However, the presence of 
even a small [si/ 2 ]^ component in the wavefunction 
would dominate the width of the decay (and lifetime). 
New detailed calculations from Grigorenko et al. using 
the three-body hyperspherical harmonics cluster model 
lowered the original constraint to Adecay < 1 keV for the 
measured lifetime [i3- Estimates from the contin¬ 
uum shell-model, assuming a sequential de cay, reported 
a similar constraint of Edecay ^ 0.5 keV [^. There¬ 
fore, it is crucial that the decay energy be further 
constrained by experimental measurements. 

In addition to the theoretical constraints on the 
decay energy, full three-body calculations of the decay 
correlations have been reported by Grigorenko et al. |17l | 
as well as by Hagino and Sagawa [s^. The results show 
that the three-body correlations are sensitive to both 
the initial wavefunction and the properties of the decay 
mechanism including the final-state nn interaction, recoil 
effects, and sub-barrier configuration mixing or rescatter¬ 
ing of the d—wave neutrons into lower i orbitals during 
the decay [sm. The predictions of the ground state 
wavefunction added additional interest since the model 
of Grigorenko et al. suggested would have a strong 
halo structure with a rms radius of the valence neutrons 
around 5.7 fm [13 and the model of Hagino and Sagawa 
showed the valence neutrons to be in a strong dineu¬ 
tron configuration . Both theoretical frameworks pre¬ 
dicted that these configurations would be manifested in 
the three-body decay correlations. One specific signature 
shown in both models is that the distribution of angles 
between the two neutrons would be peaked near 

180°. 

In this article, we extract the experimental energy and 
angular correlations in the Jacobi coordinate system from 
the three-body decay of the ground-state resonance. 
The experimental distributions are compared to a dineu¬ 
tron model, a phase space model, and the theoretical cal¬ 
culation from Hagino and Sagawa. The sensitivity of the 
experimental results to the decay mode and ground-state 
resonance energy are presented. 


II. EXPERIMENTAL DETAILS AND ANALYSIS 

Since the experimental details have already been pro¬ 
vided in Refs. [H, [s^, iHI where the ground state reso¬ 
nance and lifetime measurements were reported, only a 
brief overview is presented. An 82 MeV/u radioactive 
ion beam was produced at the National Superconducting 
Gyclotron Laboratory at Michigan State University from 
the projectile fragmentation of a 140 MeV/u ^®Ga pri¬ 
mary beam. The 2n-unbound was produced from 
a one-proton knockout of the secondary beam us¬ 


ing a 705 mg/cm^ Be reaction target. The decayed 
into -I- 2n in the reaction target. The triple coinci¬ 
dence measurement was accomplished using the MoNA- 
Sweeper setup [sol - lioji and allowed for the invariant mass 
and correlations between the decay products of to 
be determined. 

Measurements of 2n decays require the discrimination 
of the “true” 2n events from the “false” 2n background 
which is generated from a single neutron producing mul¬ 
tiple hits within the array. To remove the majority of 
false 2n events from the subsequent analysis, we applied 
the same causality cuts as used in our previous works re¬ 
porting on m, m, [ 3 ^. The causality cuts required 
the first two time ordered interactions within MoNA to 
have a spatial separation of > 25 cm and a relative ve¬ 
locity > 7 cm/ns. After applying the causality cuts, the 
false 2n component present in the experimental spectra 
is nearly negligible. 

In order to isolate the ground state decay of ^®0, only 
events with Adecay < 0.7 MeV were selected. This selec¬ 
tion was chosen to minimize any contribution from the 
decay of the first excited state of ^®0. Recent calcula¬ 
tions by Bogner et al. and Hagino and Sagawa [4l| 
indicate that the 2+ state is likely between 1 and 2 MeV. 
The Edecay specti'um with the 0.7 MeV cut is shown in 
Fig. Ha) . The full Adecay spectrum can be seen in Fig. I 
of Ref. 13 . 

The T and Y Jacobi coordinate systems, illustrated 
in Fig. [21 were used to define the angular and energy 
correlations in the three-body decay of The three- 
body correlations can be fully described by the energy. 
Ex, and angle, cos{9k), defined in each Jacobi system as, 


Ex = 


{nil + m2)kl 
2mim2 


//I \ kx * ky 

cos{9k) = 

with the Jacobi momenta kx and ky defined as. 


, m2fci — mik2 

kx = -;- 

mi -I- m 2 

^ _ m3(fci -k k^) - (mi -I- m2)ks 
^ mi + m2 -I- m3 


(la) 

(lb) 


(2a) 

(2b) 


The mass and momentum of each particle is labeled as m^ 
and ki, respectively. As depicted in Fig.|2l Ex represents 
the energy in the two-body system defined by particles 
1 and 2, while 9k represents the angle between that two- 
body system and particle 3 . Ex is often reported relative 
to the total three body decay energy (Ex). The exper¬ 
imental Ex/Et and cos{9k) distributions for the T and 
Y systems, with the causality cuts and Adecay <0.7 MeV 
criteria, are shown in panels (b), (c), (e), and (f) of Fig.jT] 


III. SIMULATIONS 

Interpretation of the data requires comparisons with 
detailed simulations of the experimental setup. A Monte 
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FIG. 1. (Color online) (a) Three-body decay energy spectrum, (b) Jacobi relative energy in T system, (c) Jacobi relative energy 
in the Y system, (d) relative velocity spectrum, (e) Jacobi angle in the T system, and (f) Jacobi angle in the Y system from 
the decay of ^®0. The experimental data are compared with simulations using three different decay modes and a half-life of 
4 ps. All results are gated on iSdecay < 0.7 MeV and have the causality cuts applied. The remaining false 2n component in the 
spectra, based on the Geant4 simulation, is shown as the solid grey area. 




T system Y system 

FIG. 2. Illustration of the T and Y Jacobi coordinate systems 
used to define the energy and angular correlations of the three 
body decay. 


Carlo simulation which included the incoming beam char¬ 
acteristics, reaction kinematics, and detector resolutions 
was utilized in the subsequent analysis (43 | . The interac¬ 
tions of the neutrons within MoNA were modeled using 
the Geant4 framework |4^ with the custom neu¬ 
tron interaction model MENATE_r . This allowed 

the results of the Monte Carlo simulation to be treated 
identically to the experimental data. 

The ground state resonance of was simulated with 
a Breit-Wigner lineshape having a resonance energy, 
and a decay width, T. The energy and angular corre¬ 
lations of the + n + n system were simulated 


assuming three different decay models: a phase-space 
decay 


[46l l47| . a dineutron decay [21|, and a decay 
model based on the theoretical calculations of Hagino 
and Sagawa [s^. The phase-space decay model simu¬ 
lates the simultaneous breakup of >-^^0-|-n -|- n as¬ 
suming that the particles do not interact during the de¬ 
cay process pll . 


e particles ac 

_ The dineutron decay model, 

described in Refs. 11 1 ^ . simulates a two-step pro¬ 
cess where a dineutron is emitted and then decays with 
an intrinsic energy defined by a nn-scattering length of 
Us = —18.7 fm [^. Lastly, a decay model was cre¬ 
ated based on the angular distribution from the full 
three-body calculations of Hagino and Sagawa (Fig. 3 of 
Ref. [11). Examples of the input cos(dfc) distributions in 
the Y Jacobi system are shown in Fig. [3] for a 5 keV reso¬ 
nance energy. As shown, the flat distribution is produced 
by the phase-space model indicating no correlations be¬ 
tween the neutrons, whereas the dineutron and Hagino 
models produces distributions that are strongly peaked 
at —1 and 1 indicating the neutrons are emitted with 
opening angles near 0° and 180°, respectively. It is im¬ 
portant to clarify that the dineutron decay model simu¬ 
lates the emission of a dineutron and does not necessitate 
the presence of a dineutron in the ground-state wavefunc- 
tion of In comparison, the correlations from Hagino 
and Sagawa are derived from a full three-body decay cal¬ 
culations which contain a dineutron configuration in the 
ground-state structure of ^®0. 
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FIG. 3. (Color online) Input cos{0k) distributions from the 
Y Jacobi coordinate system for the Monte Carlo simula¬ 
tions. The distributions from the phase-space, dineutron, and 
Hagino decay models are presented. 

The effects of a finite lifetime for were also included 
in the Monte Carlo simulation. As described in Ref. , 
a finite lifetime of in the range of picoseconds, will 
result in decreased velocity distributions of the two neu¬ 
trons emitted after the has traversed through part 
of the target. The decreased neutron velocity is repre¬ 
sented in the relative velocity spectrum, 14ei, defined as 
the difference in velocity between the neutron and the 
fragment [Fig. [TKd)]. 

IV. RESULTS AND DISCUSSION 

The results from the simulations with the three differ¬ 
ent decay modes were simultaneously fit to the experi¬ 
mental decay energy, relative velocity, and Jacobi plots 
in Fig. [TJ The ground-state resonance energy was a free 
parameter varying from 0 keV to 100 keV in the fit. The 
width of the decay, F, was fixed at 1 keV since the width 
of the resonance is completely dominated by the experi¬ 
mental resolutions [2^, [^. The influence of the life¬ 
time was also investigated with ti /2 = 0 ps and 4 ps. The 
best fit (shown in Fig. [J) of the phase-space, dineutron, 
and Hagino decay modes corresponds to Er = 15 keV, 
15 keV, and 10 keV, respectively, with ti ^2 = 4 ps. 

The experimental results were well reproduced by all 
of the models indicating that the angular-energy correla¬ 
tions are relatively insensitive to the decay mode. This is 
unexpected as previous works have shown E^/Et from 
the T system and cos(0fc) from the Y sy stem to be par¬ 
ticularly sensitive to the decay mode . For ex¬ 

ample, the input cos{0k) distribution should be strongly 
peaked at —1 for the dineutron decay, relatively flat for 
the phase-space decay, and peaked near 1 for the Hagino 
decay, as shown in Fig. [31 However, all simulations show 
the same correlation signatures when passed through the 
experimental filter. 

The insensitivity of the results to the decay mode is 
due to the uncertainty in the momentum of the at 



FIG. 4. (Color online) cos{6k) distributions from the Y Ja¬ 
cobi system for simulations using the Hagino decay mode with 
varying Be target thicknesses and decay energies for The 
results are shown (a) with and (b) without the inclusion of 
the experimental incoming beam characteristics in the simu¬ 
lations. 


the reaction point in the target and the near threshold 
energy of the ground-state resonance. The Jacobi 
variables, which describe the three-body correlations, are 
directly related to the relative momenta of the particles 
in the decay. When the decay energy is very low, the 
relative momentum between the particles is then very 
small. Thus, a small uncertainty in the momentum of 
the fragment can cause a large change in the relative 
momentum between the neutrons and fragment. This 
produces a false enhancement in the correlation between 
the two neutrons relative to the fragment which is 
manifested as the observed dineutron signatures shown in 
panels (b) and (f) of Fig.jTJ This effect is diminished with 
increasing decay energy. For example, in Fig. SKa) the 
the cos(0fc) distribution is shown for a resonance energy 
of 200 keV with the Hagino decay model (black dash-dot 
line). With the 200 keV resonance energy, the angular 
distribution peaks at cos(0fc) = 1 reproducing the overall 
feature of the input distribution from Hagino. 

The uncertainty in the momentum of the stems 
mainly from two factors: (1) the uncertainty in the po¬ 
sition at which the ^^F(—p) reaction occurs within the 
target and (2) the size and angular spread of the incom¬ 
ing beam. While the momentum, or Hp, of the outgoing 
is defined by the path of the fragment through the 
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Sweeper magnet, the momentum at the reaction point 
requires the location of the reaction in the target to be 
known so that the energy-loss within the target can be 
accurately calculated. Since the reaction location is un¬ 
known, it was assumed that the reaction took place at 
the mid-point of the target in the present analysis. In 
Fig-IHa), the effect of the target thickness is shown. For 
a 200 keV resonance, decreasing the target thickness from 
705 mg/cm^ to 0.5 mg/cm^ improves the agreement be¬ 
tween the results of the simulation and the input Hagino 
distribution. However, decreasing the target thickness 
does not improve the results when Er = 5 keV and the 
dineutron signature still prevails. This is due to the char¬ 
acteristics of the incoming beam which can produce 
uncertainties in the Bp analysis of the through the 
Sweeper magnet. Again, if the decay energy of is 
near threshold any uncertainty in the momentum is 
magnified in the three-body correlations. From the ex¬ 
periment, the beam spot on the target was estimated to 
be 5 mm x 5 mm and the incoming angular distribution 
of the beam was approximated as a Gaussian distribution 
with a mean of 0° and a = 1 mrad (2 mrad) in the dis¬ 
persive (non-dispersive) direction. These deviations from 
an ideal beam could be corrected event-by-event using 
tracking detectors upstream of the target, however such 
detectors were not available for the present experiment. 

Fig-IHb) shows the expected angular distributions as¬ 
suming ideal beam characteristics. While the simulations 
with the 200 keV resonance energy are not significantly 
altered with the ideal beam characteristics, the dineu¬ 
tron signature is no longer present in the simulation with 
a 5 keV resonance energy. A relatively flat cos(0fe) distri¬ 
bution is observed with a 125 mg/cm^ target and a slight 
peaking towards cos(0fc) = 1 is shown with a 0.5 mg/cm^ 
target with Er = 5 keV. The results from the simula¬ 
tion, therefore, indicate that an experimental measure¬ 
ment sensitive to the angular and energy correlations 
in the three body decay of will be extremely diffi¬ 
cult assuming the ground state resonance is near thresh¬ 
old. Even with an increased ^^F beam rate, the use of 
a 0.5 mg/cm^ target would pose a significant challenge 
in obtaining the required statistics to identify the decay 
mode of ^®0. While future experiments will likely en¬ 
counter a similar situation, it is important to emphasis 
that these results and sensitivity studies are specific to 
the present experiment and MoNA-Sweeper setup. 


V. CONSTRAINTS ON THE 2®0 
GROUND-STATE RESONANCE ENERGY 

Although not sensitive to the decay mode, the mea¬ 
sured three-body correlations are sensitive to the 
decay energy and, therefore, can allow for improved con¬ 
straints to be extracted through a self-consistent fit of all 
the spectra presented in Fig.[T] As mentioned previously, 
the width of the resonance was kept constant, F = 1 keV. 
The htting procedure was completed for ti /2 =0 ps and 



(keV) (keV) 


FIG. 5. (Color online) Reduced chi-squared as a function of 
^®0 gronnd-state decay energy for (a) ti /2 =0 ps and (b) 4 ps 
with the three different decay modes used in the simulation. 


TABLE I. Icr limits on the ^®0 ground-state resonance en¬ 
ergy extracted from the simultaneous fitting of the three-body 
correlations and decay energy spectrum with different decay 
modes in the simulation. 



la limit 

Model 

tl/2 = 0 ps 

ti /2 = 4 ps 

Hagino 

< 15 keV 

< 31 keV 

Phase-Space 

< 33 keV 

< 53 keV 

Dineutron 

6-42 keV 

6-53 keV 


4 ps to examine the half-life dependence of the results. 
Since the Vrei distribution [Fi^Hd)] has been shown to 
be best fit with ti /2 4 ps [32|, it was not included in 
the fitting procedure with ti /2 = 0 ps. 

The results of the simultaneous fitting procedure for 
ti /2 = 0 ps and 4 ps are presented in Fig. [5] where the 
reduced chi-squared value iXred) from each fit is shown 
as a function of the resonance energy for the three 
simulated decay modes. In all cases, it is clear that the in¬ 
clusion of the three-body correlations into the fit greatly 
enhances the sensitivity of the results to the decay en¬ 
ergy relative to our previous work where a constraint of 
TSdecay < 200 keV was extracted using only the decay 
energy spectrum . The strong correlations present in 
the input distributions of the Hagino and dineutron de¬ 
cay models (Fig. [3]) produce a strong sensitivity of the 
results at high and low decay energies, respectively. At 
low decay energies, the correlation between the emitted 
neutrons are severely overestimated by the dineutron de¬ 
cay resulting in the large Xred- Similarly, as the decay 
energy increases the correlations from the Hagino decay, 
with cos(0fe) peaking near I, become more prominent pro¬ 
ducing a rapid increase in Xred- "Fhe phase-space decay 
is the least sensitive to the decay energy since it does not 
exhibit such strong features in the correlations. 

For all three models, the Xred with the decay en- 
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ergy increasing beyond '^20 keV due to the poor fitting 
of the three-body correlations. From the chi-squared fits, 
the Ifj limits for the ground-state resonance energy 
were extracted and are shown in Table HI The results 
strongly support the conclusion that the ground-state 
resonance energy is near threshold, with upper-limits 
ranging from 15 keV to 50 keV depending on the decay 
mode and half-life. The extracted Icr limits are in excel¬ 
lent agreement with the 40 keV upper-limit reported by 
the R3B collaboration [2^. It is worth noting that the 
results using the correlations from the full three-body 
calculations of Hagino and Sagawa [s^ (which also agree 
with the detailed calculations of Grigorenko et al. 17|) 
should provide the most realistic model for the decay and 
produces Icr upper-limits of 15 keV and 31 keV for the 
cases of ti /2 = 0 ps and 4 ps, respectively. 


VI. CONCLUSIONS 

The three-body energy and angular correlations in the 
ground-state resonance decay of + n + n 

were experimentally measured using the MoNA-Sweeper 
setup. The experimental results were compared to Monte 
Carlo simulations using three different decay models: a 
phase-space model, dineutron model, and three-body de¬ 
cay model based on the theoretical calculations of Hagino 
and Sagawa Q. The experimental three-body correla¬ 
tions were well reproduced by all three of the decay mod¬ 
els indicating an insensitivity of the experimental data 
to the decay mode. This was shown to be due to the 
low decay energy of which, therefore, requires very 
precise measurements of the relative momentum of the 


and two neutrons to reconstruct the correlations. 
Monte Carlo simulations showed that the target thick¬ 
ness and tracking of the incoming beam characteristics 
largely define the experimental resolutions for measuring 
the three-body correlations. Even with improving these 
aspects of the experiment, the simulation results indicate 
that future attempts to measure the three-body correla¬ 
tions from the decay of will face a difficult challenge 
due to the near threshold ground-state resonance energy. 

Through simultaneous and self-consistent fitting of the 
decay energy spectrum and the Jacobi three-body corre¬ 
lation variables strict constraints on the decay en¬ 
ergy were extracted. While the results were dependent 
on the decay mode and half-life of a maximum Icr 
upper-limit of 53 keV was obtained. Furthermore, if one 
assumes that the correlations from the full three-body 
decay calculations of Hagino and Sagawa are cor¬ 
rect (noting that similar correlations are also predicted 
by Grigorenko et al. [13): then the Icr upper-limit for the 
ground-state is 31 keV. Additional measurements, re¬ 
quiring increased statistics, are strongly desired to pro¬ 
vide a precise measurement ground-state resonance en¬ 
ergy for 
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